The Globular Cluster Mass Function as a Remnant of Violent 

Birth 

Bruce G. Elmegreen 

IBM T. J. Watson Research Center, 1101 Kitchawan Road, Yorktown Heights, New York 

10598 USA 

bge@us . ibm. com 
ABSTRACT 

The log-normal shape of the mass function for metal-poor halo globular clus- 
ters is proposed to result from an initial power law modified rapidly by 
evaporation, collisions with clouds, and mutual cluster interactions in the dense 
environment of a redshift z ~ 5 — 15 disk galaxy. Galaxy interactions subse- 
quently spray these clusters into the galaxy group environment, where they fall 
into other growing galaxies and populate their halos. Clusters forming later in 
z ^ 2 — 5 galaxies, and those formed during major mergers, produce metal-rich 
globulars. Monte Carlo models of evolving cluster populations demonstrate the 
early formation of a log-normal mass function for typical conditions in high- 
redshift galaxies. 

Subject headings: globular clusters: general — galaxies: formation — galaxies: 
starburst — stars: formation 



1. Introduction 



The globular cluster mass function (GCMF) in prese nt-day galaxy halos is approxi 



mately log-normal with a peak Mp ~ 10 Mq (see reviews in lMcLaughlinll2003uBrodie &: Strader 



20061 ). The origin of this peaked distribution is not understood. Wherever massive dense 
clusters like globular clusters (GCs) are formed today, they have a power-law mass function 
like dN/dM cx M~^, or a Schec hter functio n with a similar power law at low mass and a 
cutoff at high mass (see review in lGielesll2009l ). Consequently, one theory for halo GCs is that 
they begin with a power-law mass function and then lose their low-raass members through 
dispersal over a Hubble ti me (iFall fc Reeslll977l : l0kazaki &: Tosalll995l : lElmegreen fc Efremov 
19971 : iFall &: Zhangll200l(). The dispersal rate wor ks out about right if the disruption process 
is thermal evaporation ( {McLaughlin fc Fall 2008 ). 



-2- 



The proble m with this model is that the GCMF do es not vary with radius in several 
nearby galaxies (ITamura et al.l l2006l : I Jordan et al.l 120071 ). and evaporation is expected to 
occur faster in the inner regions where tidal forces are larger, thereby shifting the peak 
toward highe r masses there. A rn odel in which GCs have a small range of pericenters can fix 
this pro blem (IFall fc Zhang||200in . but the resulting GC velocities disagree with observations 
in M87 JVesperini et al.lboosh . 



Wh at is expected to drive the ra dial gradient in evaporation rate is a gradient in the tidal 
density (iGieles fc Baumgardtl l2008l ). which is the average GC density inside the GC tidal 
radius. This tidal density is difficult to observe directly and is not necessarily proportional to 
the average density inside the half-light radius, which is observe d directly. The half-light den- 
sity d oes n ot correlate well with rad ial distance. Noting this, IChandar. Fall fc McLaughlin 
fl2007h and lMcLaughhn fc Fall! J2008I ) fit the GCMFs in M104 and the Milky Way to evolved 
Schechter functions for three bins of half-light density using GC evaporation rates propor- 
tional to the square roots of these densities. The results are consistent with faster evapora- 
tion, i.e., higher pe a k mas s, at greater half-light density, regardless of galactocentric radius. 
McLaughhn k Fall! J2008[ l also fit the Milky Way GCMF for three bins of tidal density using 
evaporation rates given by that density. 

An alternative model is that GCs are born with a peaked mass funct ion and then 
keep it over a Hubble time ( IVesperini et al.ll2003l : iParmentier fc Gilmorell2005l ). Log-normal 
mass functions evolve somewhat self-similarly duri ng evaporation, and their peak mass and 
width may even converge to the observed values (IVesperinil Il998l ) . The peak mass would 
also be uniform with galactocentric radius after a while. Initially peaked GCMFs could 
result from a lack of low-mass clouds in the early GC environment (IParmentier fc Gilmore 
20071 ). Another model considers variable star formation efficienci es with a greater prob- 



ability for lower cluster masses to disp erse when the gas leaves (IParmentier et al.l 12008 



Baumgardt. Kroupa fc Parmentied 120081 ). A third model is that low mass clusters were born 
with lower central conc entrations and so evaporated more quickly than high mass clusters 
( IVesperini fc Zepj |2003[ ) . There is no direct evidence for any of these models because all 
known clusters today are born with power-law mass functions. For one of these models to 
be viable, the qualitative nature of cluster formation would have to be different in the early 
Universe. 

We consider here a model where halo GCs form with power-law mass functions that are 
quickly converted into peaked functions in very den se cloudy enviro nments. The GCMF stays 
peaked and insensitive to environment thereafter (jVesperinil Il998l ) . This model involves the 
same physical process of cluster formation that is present today, i.e., gravitational collapse 
in giant gas complexes, and the cluster-forming cloud-cores are probably similar as well. 
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considering that GC densities, masses, and IMFs are not unusual. What differed in young 
galaxies was a much dens er and more turbulent interstellar medium (ISM) than we have i n 
main galaxy disks today ( iForster Schreiber et al.ll2006l : iGenzel et al.ll2006t iLaw et al.l 120071 ). 
GC formation in z > 10 dwarf galaxies was discussed by iBromm J2004|)T bulge-GC form ation 
in 2 ~ 2 galaxy clumps was discussed by IShapiro. Genzel fc Forster Schreiberl (12010 ), an d 
simulations of GC formation in young galaxy disks were made by Kravtsov fc GnedinI ( 2005 ). 
The main point here is that when the density is high and the motions are fast, clusters should 
frequently collide with cloud clumps and other clusters during their first several hundred 
Myrs. These collisions destroy the lowest mass clusters and produce a log-normal GCMF. 



2. Collisional Dispersal of Clusters in High Density Environments 



Clusters heat up or disperse completely when their potential energy changes quickly. 
This can result from rapid gas expulsion or from collisions with dense clouds and other 
clusters. In a star-forming environment, there are many opportunities for collisions. The 
cluster dispersal rate therefore starts high after birth, and it decreases as the star-forming 
cloud dissipates and all the clusters and dense cores drift away . We have studied time- 
dependent cluster collisions elsewhere ( lElmegreen fc Hunterll2010l ). Here we consider a time- 
average collision rate and ask what it has to be in order to convert a power-law initial cluster 
mass function into a peaked mass function during the starburst phase of a young galaxy. 



The collisional disruption rate is taken from iGieles et al.l (120061 ) . who use a cluster 
evolution equation dM/dt = — M/tdis with a disruption time 



tdis = 37- Myr; 

^nPn 



En is the average column density of a collision partner, in Mq pc ^, pn is the average density 
of collision partners in the neighborhood, in Mq pc~^, and M is the cluster mass. Gieles et 
al. use En = I7OM0 pc'^ for typical GMCs and Pn = O.OSMq pc'^ for the ambient ISM. 
Then E^p^ ~ 5.1M2 pc"^ and tdis ~ 7.3M°-^^ Myr, which is 8.2 Gyr for M = 10^ Mq. 

The disruption time is less at higher density. If we require disruption of M < 10^ Mq 
clusters within the time span of t he gas-rich phase of a young galaxy disk, which may be 
500 Myr (e.g., iTacconi et aPboOsI ). then tdis = 500 Myr and SnPn ^ 80M^ pc"^ - a factor 



of 16 larger than the local ISM value used by Gieles et al.. This factor may be accounted 
for by a higher average density in the cluster environment and a higher average column 
density for likely collision partners. For example, pn c ould be ~ 10 times higher in a galaxy 



with 50% of the disk mass in the form of clumpy gas (ITacconi et al.ll2010l ). and En could be 



-4- 



10^ Mq pc ^ or more for massive molecular cores and other clusters. If we require that 



^dis equals 10 times the dynamical time of the ISM, which is 10/ {Gpnf'^, then Snp°-^ has to 



equal 280 Mq pc~ for M = 10 Mq. Even this is reasonable for an extremely gas-rich 
and dense disk where might be ~ 1 Mq pc~^ (~ 30 atoms cm~^) in regions one kpc in 
size. The clumpy gala xies observed at z ~ 2 are a good example of such an environment 
( lElmegreen et al.l 120091 ). The ISM density should be even larger at z ~ 5 — 15, where many 
halo GCs are likely to have formed. Young galaxies are denser than today's galaxies because 
the Universe was denser, by the factor (1 + 2;)^ ^ 1000 for the redsh ift of interest here. Young 
galaxies are also smaller than today's galaxies (lOesch et al.ll2010[ ) and either represent the 
inner (dense) parts of today's galaxies or separate objects that merged into dense spheroids 
over time. 



2.1. Monte Carlo Simulations with Constant Cluster Birthrate 



The time evolution of cluster populations was determined for Monte Carlo simulations 
with constant cluster birthrates and cluster disruption rates given by 



dM/dt = -M/Uis ; Uis = ^M^; 



(2) 



^ = 37/ (SnPn) could be ~ 0.01 or smal l er for units of column density in Mq pc ^ and 
density in Mq pc~^. In the iGieles et al. I J2006h model, 7 ~ 0.6; we take 7 = 0.62 to be 



consistent with the evaporation model i r ilLarners et al.l (120051 ) . We also ran simul ations with 
7 = 1, which is appropriate for ISpitzerl (119871 ) evaporation and for ISpitzerl (119581) coUisional 
disrupti on with a weak mas s-rad ius relat i on, as observed by lBastian et al.l (120051 ) and others. 
That is, iGieles et al.l (120061 ) and ISpitzen (119581 ) both derived tdis oc pd/ (X'nPn) for coUisional 
disruption with cluster internal density pci- Bastian et al. used pci oc M°-^^, but it is also 
possible that pci oc M. Both coUisional disruption and cluster evaporation should be rapid 
in the high density environments of young GCs. CoUisi onal disruption probably dominates 
evaporation in the disk environment (IGieles et al.ll2006l ). 



Values of C, that give a mass distribution peak Mp when the cluster population has an 
age T may be determined from the relation ^MJ = T. For T = 1000 Myr and Mp = 10^ Mq, 
^ = 0.01 when 7 = 1, and ^ = 0.79 when 7 = 0.62. These ^ are reasonable for starburst 
conditions in young, dense galaxies. 

In the simulations, clusters masses were randomly chosen from an initial dn/dM cx M~^ 
mass function that extends from a minimum mass Mmin = 10 Mq to a maximum mass of 
10® Mq. The cluster formation rate is one cluster per time step, dt, measured in Myr. At 
each time step, equation ([2]) is applied to every cluster, and every cluster mass is reduced 
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accordingly. Clusters with masses dropping below O.lMmin are not followed. The results 
after 10 Gyr are shown on a plot of log M versus log T for each current cluster mass M and 
age T (Fig. [T]). The age is defined to be the current time in the simulation minus the time 
of formation of the cluster. We are interested in the distribution of cluster mass at an age of 
~ 500 Myr or so, when the starburst ends and the clusters begin to scatter. This distribution 
can be read off the plot at T = 500 Myr. The plot extends to T = 10 Gyr for comparison. 

Figure [T] shows log M — log T distributions of cluster populations for various values of 
^ and for 7 = 0.62 on the left and 7 = 1 on the right. All of the distributions have about 
the same form with more or less attrition over age in cases with low or high ^, respectively. 
The upper envelope of the cluster mass comes from the size-of-sample effect, which states 
that the most massive cluster in an dn/dM ~ distribution is directly proportional to 
the number of clusters in the sample. This maximum increases linearly with T on a plot 
like this because the time interval AT increases linearly with T for equal intervals in logT. 
The lower envelope of the distributions increase as M oc T^/''' because that is the mass at 
which tdis = T. The dashed green hne is Mmm = 10 Mq. For 7 = 1, the upper and lower 
envelopes are parallel so the mass range of the mass function at any one age remains about 
constant. For 7 = 0.62, the bottom envelope increases faster than the top envelope. Two 
cluster formation rates {1/dt) are used to give good numbers of points at different ^. 

Lower ^ makes the characteristic mass of the clusters higher for a given age. Figure [1] 
has vertical and horizontal dotted lines at fiducial markers T = 1000 Myr and M = 10^ Mq. 
At T = 1000 Myr with 7 = 0.62, ^ = 0.1 has a mass distribution centered slightly higher 
than M ~ 10^ M©, while ^ = 1 has a mass distribution centered slightly lower than this. 
The value ^ = 0.5 gives a mass at the peak of the function Mp ~ 10^ Mq, as predicted 
approximately above. For 7 = 1, a value of ^ ~ 0.01 gives this Mp at T = 1000 Myr, as also 
predicted. 

Figure [2] shows log-normal fits to the cluster mass functions versus cluster age in four 
age intervals, logT = 2 — 2.5, 2.5 — 3, 3 — 3.5, and 3.5 — 4 for T in Myr. The trend of 
increasing mean log M with increasing age (Figure [1]) is evident in Figure [2l The scatter in 
the diagram is from the small numbers of clusters in the age intervals (50 to 500, depending 
on C, and 7). The solutions give the desired Mp ~ 10^ when T ~ 10^ Myr if ^ is small. 

2.2. Monte Carlo Simulations with Cluster Birth for 500 Myr 

Figure [3] shows another case. The blue dots come from previous models (Figured]). The 
red crosses are for a model where cluster formation occurs for the first 500 Myr and cluster 
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disruption uses equation (|2]), but then star formation stops and the disruption rate decreases 
in proportion to the remaining cluster mass as 



Mtot(t = 500 



(3) 



for times between 500 Myr and 1000 Myr. This model simulates normal cluster formation 
and disruption during 500 Myr when dense gas is present, followed by no cluster formation 
and diminished disruption for another 500 Myr after the gas disperses. Only the log-normal 
GCMF remains; fits are shown in Figure HI Figure [2] shows black symbols for the peak masses 
and dispersions in this second model at T = 2.75 — 3.25 Myr; they should be compared with 
the black curves, which have the same 7 and C,- 



Discussion 



We propose that old halo GCs formed by normal processes much like clusters form today, 
but that all of these processes occurred in young galaxies at very high densities, higher than 
today's ISM densities by factors of 10 to 100. Then coUisional disruption and evaporation 
was rapid enough to produce a peaked mass function in the young galaxy. Subsequent 
evolution by evaporation in lower density environments preserved this function, although 
slight changes in peak mass and width probably occurred. 

Given this basic model, there are many possibilities for the delivery of these clusters 
into modern galaxies. The dense galaxies they formed in are most likely not the same as the 
galaxies or inner parts of the galaxies that currently host them. Their presence in galaxy 
halos today implies that they were delivered to the host in a non-dissipative way, perhaps by 
infall along with most of the host's other baryons. For example, the oldest GCs could have 
formed in small dense galaxies that form ed as condensations in larger-scale cold gas flows . 



The cold flows rn ake today's spiral disks (jPekel et al.ll2009t lAgertz. Teyssier fc Moord 12009 



Keres et al.ll2009l ). and the GCs along with other condensations in the flow populate modern 
halos. The small dense galaxies could also have collided before they entered the modern 
galaxy's potential well, freeing up the GCs which would then enter the well as free-floaters 
along with the other material. Some smal l galaxies are still adding GCs to rnodern halos 



[e.g.. iCarraro et al. 1120071 : ICao et all 120071 ICasetti-Dinescul l2009l : ISmith et al.ll2009[ ). The 
oldest GC populations in elliptical galaxies presumably formed in small dense galaxies and 
fell into spiral halos too in the same way, but then ended up in ellipticals after major mergers 
of the spirals. These oldest GCs would be the blue, metal-poor populations in spiral and 
elliptical halos. GCs fo rmed during major mergers would be redder and more metal-rich 
JSrodie fc Stradeilbood ). 
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The uniform properties of old halo globular clusters (color, density, IMF, peak GCMF 
mass, metallicity) follow from this model if we postulate that GCs are the first examples of 
star formation at high enough metallicity to make a normal IMF. Then stellar evolution will 
not remove excessive amounts of gas and cause the cluster to come unbound. Top-heavy 



IMFs, such as those thou ght to produce the carbon-enhanced metal-poor stars ( iTumlinson 



20071 : iKomiya et al.l 120091 ). lose too high a fraction of their mass during stellar evolution to 
keep a cluster bound. Boundedness would require them to occupy the nuclei of small galaxies 
so they can retain their stellar wind material in the galactic potential well. Then contam 
ination from subsequent ge nerations of stars would also enrich them (e.g., iMarcolini et al. 



20071 : iBailin fc HarrisI 120091 ). For clusters born in small dense disks, the uniformity of both 
the clusters and the GCMFs suggest that the first epoch of normal IMFs occurred when 
galaxies were still small and dense, and this was long before today's spirals were assembled. 

The bulge GCs of modern spirals could have a similar origin as the halo GCs, but there 
is no similar constraint on the constancy of their GCMFs over a range of galactocentric 
radii. Also, the bulge GCs formed in metal-rich environments and they are now located 
deep in the potential wells of their current hosts. Thus they need not have formed in 
separate small galaxies. They could have for med in their hosts and settled to the cente r 
after energy dissipation. This is the model by IShapiro. Genzel &: Forster Schreiberl ( 120101 ). 
They proposed that bulge GCs formed in the dense clumps of young massive z ~ 2 galaxies, 
much like we propose that metal-poor GCs formed in the dense clumps of young and small 
2; ~ 5 — 15 galaxies. The local environment is about the same for each, namely high-density 
clumpy gas, so collisional disruption would have occurred quickly for the Shapiro et al. model 
too. The densities in the massive clumps of 2; ~ 2 spirals are not expected to have been 
as high as the densities in the disk clumps of smaller galaxies at higher redshift, however. 
Thus collisional disruption of low mass GCs may not have been as rapid for the Shapiro et 
al. model. 

Our model accounts for the number of halo GCs today if a high fraction of early star 
formation made bound clust ers. Such a high fraction is expected when the ISM pressure and 
velocity dispersion are large (lElmegreenll2008[) . The space density of GCs today is ~ 8 Mpc~^ 
( iPortegies Zwart fc McMillan! 120001). There could have been a factor of ~ 2 more when they 
formed, considering evaporation (jVesperinil Il998[ ) . but blue halo GCs are only half of the 
total. The co-moving s pace density of clumpy galaxies at 2; ~ 1 — 2 is rig ~ 2 x 10^^ Mpc~'^ 
( lElmegreen et al.ll2007l ). Each galaxy has ~ 5 clumps of stellar mass ~ 10^ M©, which would 
be iVc ~ 5 X 10^ clusters of mass ~ 10^ Mq each. Thus the space density of m assive clusters 
would have been n„Nc ~ 10 Mpc~^, with a factor-of-3 uncertainty either way. IShapiro et al. 



( 120101 ) got the same result a different way. If metal-poor GCs formed in younger, smaller 

gal 



galaxies, then the M ? cosmological galaxy mass function means that each log interval of 
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galaxy mass contains the same total mass. Thus we would get these same ~ 10 GC per 
Mpc^ for GC formation in galaxies 10 or 100 times less massive at higher redshifts. 

Modern examples of this model may occur in dense galactic nuclear regions where we 
predict a rapid destruction of low mass clusters because of heightened collision and evapo- 
ration rates. 

Acknowledgements: I am grateful to Enrico Vesperini and Tanuka Chattopadhyay for 
discussions about the GCMF, and to Dean McLaughlin for comments on the manuscript. 
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Fig. 1. — Distribution of cluster masses M and ages T after 10^ Myr of evolution with 
continuous cluster formation in an M~'^dM mass function and continuous disruption or 
evaporation according to equation (2) with values of 7 and ^ indicated. The star formation 
rate is 1/dt for model time step dt in Myr. The horizontal dashed line is the lower limit 
to the mass of a formed cluster; clusters get lower mass over time because of the assumed 
disruption. The blue dotted lines indicate fiducial markers for age and mass where we expect 
a GC population should be after the starburst phase of a young galaxy. The clusters are 
presumed to be scattered after ~ 1 Gyr and then they evolve mostly by slow evaporation in 
the halos and bulges of these and other galaxies. 
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Fig. 2. — Fits for peak masses Mp and dispersions a of log-normal cluster mass functions 
for 0.5 log-age intervals. Different curves are for different ^, which is the coefficient in 
equation (2) for the dispersal time. The peak mass increases with age because of preferential 
dispersal of low mass clusters. Strong cluster disruption, corresponding to small values of ^, 
can produce a peaked GCMF with Mp ~ 10^ Mq after ~ 1 Gyr, as indicated by the dotted 
lines. Diamonds give Mp and a values for the models shown as red crosses in Fig. 3. 
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Fig. 3. — Masses and ages of cluster populations after 10^ Myr with continuous cluster 
formation and dispersal (blue dots, as in Figure 1). Surviving clusters (red crosses) for a 
second model in which cluster formation and full-rate disruption end after 500 Myr, and 
then cluster disruption diminishes in proportion to the total mass of the remaining clusters. 
This second model is viewed after 10^ Myr. It represents a case where clusters are formed 
and disrupted during a Gyr time span in a small dense galaxy at high redshift. The clusters 
end this phase with a log-normal mass function and presumably scatter into other forming 
galaxies over time, preserving this mass function. 
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log M log M 

Fig. 4. — Mass distribution functions of the red crosses shown in Figure 3, along with log- 
normal fits. The 7 = 0.62 case has a tail of low mass clusters for the assumed epoch (1 Gyr) 
and disruption rate (^), but both are reasonably peaked at around Mp ~ 10^ Mq, making 
this process viable as a mechanism to inject GCs with this "initial" mass function into the 
converging flows of the early universe. 



